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Symmetry of spin excitation spectra in the tetragonal paramagnetic 
and superconducting phases of 122-ferropnictides 
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We study the symmetry of spin excitation spectra in 122-ferropnictide superconductors by compar- 
ing the results of first-principles calculations with inelastic neutron scattering (INS) measurements on 
BaFci 85C00.15AS2 and BaFe^ 9iNio.o9As2 samples that exhibit neither static magnetic phases nor structural 
phase transitions. In both the normal and superconducting (SC) states, the spectrum lacks the three- 
dimensional (3D) 42/ m screw symmetry around the (\\L) axis that is implied by the lA/mmm space group. 
This is manifest both in the in-plane anisotropy of the normal- and SC-state spin dynamics and in the out-of- 
plane dispersion of the spin-resonance mode. We show that this effect originates from the higher symmetry 
of the magnetic Fe-sublattice with respect to the crystal itself, hence the INS signal inherits the symmetry of 
the unfolded Brillouin zone (BZ) of the Fe-sublattice. The in-plane anisotropy is temperature-independent 
and can be qualitatively reproduced in normal-state density-functional-theory calculations without invoking 
a symmetry-broken ("nematic") ground state that was previously proposed as an explanation for this effect. 
Below the SC transition, the energy of the magnetic resonant mode co^^^, as well as its intensity and the SC 
spin gap inherit the normal-state intensity modulation along the out-of-plane direction L with a period twice 
larger than expected from the body-centered-tetragonal BZ symmetry. The amplitude of this modulation 
decreases at higher doping, providing an analogy to the splitting between even and odd resonant modes in 
bilayer cuprates. Combining our and previous data, we show that at odd L a universal linear relationship 
^cOj-es ^ 4.3 k^T^ holds for all the studied Fe-based superconductors, independent of their carrier type. Its 
validity down to the lowest doping levels is consistent with weaker electron correlations in ferropnictides as 
compared to the underdoped cuprates. 
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I. INTRODUCTION 
A. Electronic symmetry-broken states 

The ground state of a paramagnetic metal naturally in- 
herits all symmetries of its underlying crystal structure. This 
generally applies to the single-particle Bloch states in the 
periodic potential of the lattice and can be generalized to 
the spectra of particle-hole excitations or collective modes, 
such as phonons. However, several mechanisms may lead to 
spontaneous breaking of this crystal symmetry as the system 
is driven by a change of some control parameter (e.g. tem- 
perature, electron doping, or pressure) towards an ordered 
ground state. For such symmetry breaking to occur, both 
electron and lattice degrees of freedom are often required,^ 
as in magneto-structural or charge-density-w^ave^""^ transi- 
tions. Occasionally, though, the electron degrees of freedom 
alone are sufficient to lead to an instability, w^hile the lat- 
tice only adjusts itself to the new^ ground state, offering 
little contribution to the overall energy gain.^ The most 
prominent examples of such electron-driven instabilities are 
spin-density-w^ave (SDW) transitions,"^'^'^ at w^hich a mag- 
netic ordering w^avevector is spontaneously chosen out of 
several equivalent Fermi surface (FS) nesting vectors, or 
Pomeranchuk instabilities that spontaneously lov\^er the FS 
symmetry. '^"^^ 

In the special case of so-called ''electronic nematic'' 



phases, ^^~^^ only the rotational symmetry of the electron 
subsystem is reduced, whereas the translational symmetry 
and, hence, the size of the Brillouin zone (BZ), are pre- 
served. Such states have been extensively studied in quasi- 
two-dimensional compounds, such as Sr3Ru207 (Ref. 17) or 
underdoped YBa2Cu3 05+^.^^"^^ Recently, electronic nematic 
phases have been also suggested for various iron-arsenide 
superconductors^^"^^ (for the latest reviews, see Refs. 36 
and 37). 

In the following, we present the results of first- 
principles calculations and inelastic neutron scattering 
(INS) measurements of the spin excitation spectra in 
the normal and superconducting (SC) states of slightly 
underdoped BaFe2_jfNi^As2 (BFNA) and optimally doped 
BaFe2_;cCo^As2 (BFCA) single crystals, which belong to the 
so-called 122-family^^"'^^ of ferropnictides. "^^""^^ Although 
the crystal structure of both compounds retains its body- 
centered-tetragonal (bet) I4/mmm symmetry (Fig. 1) down 
to the lowest measurable temperatures, the low-energy spin 
response both in the normal and SC states has a lower 
symmetry in the reciprocal space,^^"^^ as sketched in Fig. 2, 
which corresponds to the unfolded BZ of the Fe-sublattice 
(Fig. 3). This unfolded zone is often introduced to simplify 
the band-structure description of the iron pnictides,^^'^^ 
but is usually considered only as a theoretical abstraction, 
because any realistic band structure of these systems is 



6 42/m 



MImmm 




Fig. 1 (color online). The reciprocal-space structure of the 
I4/mmm crystal. The BZ polyhedron of BaFe2As2 is drawn at 
the left in solid black lines, and two more such polyhedra are 
drawn to illustrate the 3D stacking of the Brillouin zones. Two FX 
vectors are shown by dashed arrows: The SDW vector of the parent 
compound Q^fm = (22-'^) ^^^ ^^^ in-plane projection Q|| = (||0). 
Symmetry axes are denoted by dash-dotted lines. 



certainly affected by the pnictogen atoms that low^er the 
symmetry of the direct lattice and, consequently, introduce 
an additional translational symmetry in the reciprocal space 
due to the BZ folding. Nevertheless, as we will demonstrate 
in the following, the absence of any appreciable magnetic 
moment on the pnictogen atoms allows for a much simpler 
description of the dynamical spin susceptibility, which expe- 
riences no structural folding and hence does not acquire the 
additional reciprocal-space symmetry expected in the back- 
folded tetragonal (structural, nonmagnetic) BZ. Therefore, 
as far as the magnetic fluctuations in the paramagnetic state 
of ferropnictides are concerned, the unfolded description of 
the spectrum becomes physically justified. 

Due to the three-dimensional (3D) character of the 122- 
systems, manifest both in their electronic structure"^^"^^ and 
in the substantial out-of-plane magnetic coupling in their un- 
doped (parent) compounds, ^"^"^^ the missing symmetry op- 
eration is essentially three-dimensional, involving all three 
crystallographic coordinates. It corresponds to the 42 /m 
screw symmetry around the (||^) axis,^^ shown in Fig. 1, 
and is equivalent to a product of a 90° in-plane rotation 
around the F point and a translation by the reciprocal lattice 
vector G = rr = (101). In the following, we will show that 
the clear absence of such screw symmetry — a conjectured 
3D analog of the electronic nematicity — can indeed be ob- 
served in the spin-excitation spectrum already in the normal 
(paramagnetic) state, both along the out-of-plane and along 
the in-plane directions of the reciprocal space. In this re- 
spect, our experimental data are in qualitative agreement 
with recent reports of anisotropic in-plane excitations seen 
both in the magnetically ordered^^'^^ and paramagnetic^^"^^ 
states. The latter were previously associated with "spin 
nematic correlations". However, a comparison with normal- 
state density-functional-theory (DFT) calculations presented 
in section II A shows good agreement between the calculated 
and measured susceptibilities, leading us to an alternative 



explanation for the lowered symmetry of the spin-excitation 
spectrum that does not require a symmetry-broken ground 
state or proximity to a quantum critical point. Instead, it 
turns out to be a direct consequence of the crystal structure 
with two Fe atoms per primitive unit cell, in which the crys- 
talline lattice that determines the BZ geometry has a lower 
symmetry than its Fe-sublattice, which is responsible for the 
magnetism. 
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B. Reciprocal space of the 122-ferropnictides 

To set the scene, in Fig. 3 we summarize some of the pos- 
sible coordinate systems and reciprocal-space notations that 
can be introduced in the 122-compounds. The figure shows 
five different Brillouin zones in the reciprocal space (right) 
and their respective primitive unit cells in direct space (left) . 
It is natural to consider two BZ types: unfolded, i.e. corre- 
sponding to the Fe-sublattice only, and folded, which takes 
full account of the remaining nonmagnetic atoms in the unit 
cell. Because of the higher symmetry of the Fe-sublattice 
with respect to the crystal itself, the unfolded zones have 
twice larger volume than their folded counterparts. Next, 
one can also distinguish between the nonmagnetic and mag- 
netically folded BZ, which correspond to the normal and 
SDW states, respectively. As a result, we end up with four 
different direct-space lattices, reciprocal-space coordinate 
systems, and BZ geometries that can be naturally intro- 
duced in the 122-compounds: (a) unfolded tetragonal (Fe^); 
(b) body-centered tetragonal (Fe2); (c) unfolded magnetic 
(Fe2); (d) doubly-folded magnetic (Fe4). The formulas in 
brackets give the number of iron atoms in the primitive unit 
cell. In addition. Fig. 3 (e) shows the simple tetragonal unit 
cell (Fe4) that defines the reciprocal-space notation used 
in the present paper, ^^ but does not represent a primitive 
unit cell of the crystal. In the following, we will concentrate 
on the normal (paramagnetic) state, and therefore will be 
mainly interested in the nonmagnetic (folded or unfolded) 
BZ. 

The 3D stacking of the I4/mmm tetragonal Brillouin 
zones with the dimensions of — x ^ x — is illustrated 
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Fig. 2. A sketch illustrating the symmetry of spin excitations: (a) as 
expected from the BZ symmetry in the absence of matrix elements; 
(b) as actually observed experimentally in doped 122-compounds. 
The surfaces schematically represent constant-intensity contours of 
the magnetic INS response. The center of each panel corresponds 
to the r point. Note that despite the lower symmetry in (b) due to 
the absence of the 42/m screw around {\\l^), the four-fold (4/m) 
rotational symmetry around (OOL) is preserved. 
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Fig. 3 (color online). Different primitive unit cells in direct space (left) that can be introduced in 122-ferropnictides and their respective 
Brillouin zones (right): (a) unfolded tetragonal BZ of the Fe-sublattice with one Fe atom per unit cell (Fe^); (b) structural body-centered- 
tetragonal BZ that corresponds to two iron atoms per primitive unit cell (Fe2); (c) unfolded magnetic BZ that corresponds to the 
magnetically ordered Fe-sublattice in the SDW state (Fe2); (d) doubly-folded magnetic BZ that results if both the lattice and magnetic 
structures are taken into account (Fe4); (e) one of the most commonly used and experimentally convenient reciprocal-space coordinate 
systems that corresponds to the BZ of a simple-tetragonal direct lattice with the parameters of the real bet crystal — the notation that we 
also adopt for the present paper. ^^ Use the checkboxes at the right to toggle the visibility of individual graphic layers. 



in Fig. 1 and is valid both for the momentum (k) and mo- 
mentum transfer (Q) spaces. In our notation, the quasi- 
tw^o-dimensional (2D) w^arped hole- and electronlike FS 
cylinders"^^'"^^'^^"^^ are centered around FAZ and XPX sym- 
metry axes along the zone boundaries, respectively. The INS 
data presented in this paper were measured in the vicinity 

of tw^o rX w^avevectors that are show^n by dashed arrow^s: 
The magnetic ordering v\^avevector of the parent compound 
Qafm — (2 2-'-) ^^^ ^^^ in-plane projection Qy = (||0). Note 
that the tw^o vectors are equivalent in a tetragonal system 

modulo the reciprocal lattice vector G = rr = (101), be- 
cause Qafm ~G = (-||0) ~ (||o). This equivalency is oblit- 
erated, however, by the magnetic order in the orthorhombic 
phase that selects Qafm ^s the preferred SDW vector. It is 
difficult to understand the out-of-plane component of this 
SDW ordering wavevector in a simple (geometric) nesting 
picture because of the equal nesting conditions at Qafm ^^^ 
Q||, imposed by the 42 /m screw symmetry. But as we will 
subsequently show in section II A, a more rigorous calcula- 
tion of the Lindhard function, taking into account the orbital 
matrix elements, is sufficient to resolve this dilemma. 

The crystal symmetry axes are shown in Fig. 1 by dash- 
dotted lines. In particular, the 42 /m screw symmetry along 
the XPX axis appears only in the bet BZ with 2 Fe atoms 
per primitive cell as a result of folding, but is found neither 
in the unfolded BZ corresponding to the Fe-sublattice be- 
cause of the missing (101) translation, nor in the magnetic 
BZ because of the spontaneously broken 4-fold rotational 
symmetry in the SDW or orthorhombic phases (see Fig. 3). 



It will be especially important for our discussion because of 
its insensitivity to electronic twinning of the crystal, i.e. the 
presence of domains with different orientations of the spon- 
taneously symmetry-broken electron states in samples with 
in-plane anisotropy or under the assumption of electronic 
nematicity. In contrast, the breaking of the 4-fold rotational 
symmetry around the FAZ axis cannot be directly observed, 
unless the sample is electronically detwinned, which can 
be achieved by the application of uniaxial pressure^^^'^^'^^ 
or an external magnetic field. "^^ It is also essential that the 
42 /m symmetry axis coincides with the Q-space location 
of the spin excitations found in INS experiments, which 
allows us to compare the magnetic intensities along this 
direction. These excitations, which constitute the subject 
of the present study, originate from the nested hole- and 
electronlike Fermi surfaces"^^"^^'^^"^^ and survive even in the 
overdoped regime,"^^ i.e. well above the onset of the static 
SDW order in the phase diagram. 

C. Normal-state spin-excitation spectrum 

The normal-state spin dynamics of 122 Fe-based super- 
conductors is dominated by an intense branch of low- 
energy spin fluctuations in the vicinity of the commensurate 
Q = (^L) wavevector. It is characteristic for a nearly an- 

tiferromagnetic (AFM) metaF^ and can be well described 
within an itinerant framework.^^'^^'^^'^^ At higher energy 
transfers, the spin excitations exhibit a dispersion that has 
an anisotropic cross-section within every L = const plane. 



This has been evidenced in time-of-flight (TOF) experiments 
covering odd, even, or half-integer L values. ^^'^^ The ob- 
served similarity to the magnetic parent compound^^^'^^'^^ 
served as a starting point for the proposed symmetry-broken 
("electronic nematic") ground state. 

Caution has to be taken, however, since in the struc- 
tural BZ (Fig. 1) the orthogonal XT and XZ vectors lying 
in the k^ky plane (which for L = correspond to the max- 
imal and minimal spin-wave velocities, respectively) are 
not equivalent. Indeed, the different shapes of the holelike 
barrels that alternate in a checkerboard manner, as seen in 
angle-resolved photoemission (ARPES) maps at a fixed exci- 
tation energy,^"^ confirm the significance of this difference. 
Moreover, electronic band structure calculations within the 
tetragonal phase yield elliptical in-plane cross-sections of 
the electronlike FS sheets around the X point,^^ which obvi- 
ously do not by themselves imply any anisotropy between 
the (110) and (110) directions, because the ellipse rotates by 
90° when shifting to the nextX point. Therefore, discussions 
of the in-plane anisotropy in 122-compounds necessarily re- 
quire consideration of the full 3D band structure, including 
the out-of-plane dispersion of the spin response along L . If 
the observed ellipticity followed the I4/mmm symmetry of 
the crystal, then the X -centered intensity pattern in the spin 
susceptibility would be rotated by 90° at odd L with respect 
to even L values because of the 42 /m screw symmetry, as 
illustrated in Fig. 2 (a) . On the contrary, the absence of this 
symmetry in the spin-excitation spectrum may lead to the 
same orientation of the ellipse at all L and to the doubling 
of the period of intensity modulation along (||l), as shown 
in Fig. 2(b). 

In order to discriminate between these two possibilities, 
we performed triple-axis INS measurements in the (^HK [H+ 
K]) scattering plane, thus avoiding the L -integration that 
is pertinent to the TOF method. A direct comparison of 
the transverse and longitudinal scans around the (||l) and 
(— 1|0) wavevectors, presented in section IIIC, shows that 
the excitation spectrum indeed does not fully follow the 
crystal symmetry, but inherits it only from the magnetically 
active Fe-sublattice. This consequence of the material's 
crystallography per se does not imply any spontaneously 
symmetry-broken states in direct space. Moreover, the van- 
ishing L -dependence of the anisotropy ratio indicates that 
the structural contribution to the ellipticity (originating from 
the folded FS geometry) is not detectable within our experi- 
mental accuracy. 



D. Superconducting spin-resonance mode 

The magnetic resonant mode is the most prominent 
signature of superconductivity in the spin-excitation spec- 
trum of several unconventional superconductors, such as 
single -layer "^^ and bi-layer cuprates^^"^^ or heavy-fermion 
systems. ^^ In this respect, the Fe-based superconduc- 
tors are no exception: A resonance was found in hole- 
doped Bai_^K^Fe2As2,^'^ optimally electron-doped BFCA 
and BFNA,^^'^^ underdoped BFCA,^^ iron chalcogenides 
FeTei_^Se^,^^'^^~^^ and more recently in polycrystalline 
LaFeAsOi_^F^.^^'^^ 



The resonant mode carries information about the symme- 
try of the SC gap — d-wave in cuprates^^ and s±-wave in the 
iron arsenides. ^^'^^^'^^^ The resonance in Fe-based supercon- 
ductors shares various common aspects with cuprates, such 
as its abrupt intensity evolution below T^, and the fact that 
it is always observed at an energy cOj-es below the particle- 
hole continuum that sets in at twice the SC gap a.^^^'^^^^ 
However, there are also differences: In BFCA, the tempera- 
ture evolution of co^es is BCS-gap-like, and no signature of a 
pseudogap has been found. ^^ 

In section IV B, we will compare two further aspects of 
the resonant features in both systems. First, due to the 
intra-bilayer coupling, bilayer cuprates exhibit two resonant 
modes characterized by odd and even symmetries with re- 
spect to the exchange of CUO2 layers within a bilayer unit, 
as reported for the YBa2Cu305+^ and Bi2Sr2CaCu208+5 
families. ^^^~^^^ These modes show intensity modulations 
with L, antiphase with respect to each other, as well as 
different but L -independent resonance energies. Although 
distinct resonance energies for even and odd L were ob- 
served in BFNA^^'^^^ (see also section IIIB), a comparison 
to the cuprates has not yet been drawn, because due to the 
equally-spaced FeAs layers, two distinct resonant modes are 
not expected. Second, a linear relationship between cOj-es 
and Tc has been extensively discussed for cuprates, and a 
ratio of co^^^/k^T^ ^ 5.3 has been established for the odd 
resonance, for doping levels not too far from optimal. ^^"^ 
However, progressive deviations have been noted with 
underdoping,^^^'^^^ a violation was reported for single-layer 
HgBa2Cu04+3,,^°^ and there is an ongoing controversy about 
the situation in electron-doped cuprates. ^^^'^^^ In contrast 
to this, as we will show in section III D, a similar linear rela- 
tionship cOj-es/^B^c ^ 4.3 is universal among all the studied 
Fe-based superconductors, over the entire phase diagram 
and independent of their carrier type, and holds down to the 
lowest doping levels. This means that the coupling strength 
(as opposed to co^es) very weakly depends on doping. 



II. RESULTS OF FIRST-PRINCIPLES CALCULATIONS 

A. Normal-state Lindhard function 

We start by presenting DFT calculations of the Lindhard 
function^ ^^"^^^ 

(O ^= ^ V /n-(k + Q)-A(k) 



j^jE,.(k+Q)-F„(k)- 

where ^„(k) is the energy of the n-th band, |k, n) is the cor- 
responding wave function, /„(k) is the Fermi function, and 
a± are Pauli matrices. These calculations were performed 
starting from the tetragonal non-magnetic state for the ex- 
perimentally determined atomic positions. ^^'^ The chemical 
doping was included in the virtual crystal approximation. 
Further details of the calculations can be found in Ref. 115. 
The surface plots of the static susceptibility ;^o(Q^ ^ ~^ 
0) in the undoped BaFe2As2, 10% Co-doped (electron- 
overdoped), and 40% K-doped (optimally hole-doped) com- 
pounds are shown in Fig. 4 for L = and L = 1 together 
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Fig. 4 (color online). The Lindhard function /o(Qj ^)^ resulting from DFT calculations in the undoped (top), 10% Co-doped (electron- 
overdoped, middle), and 40% K-doped (optimally hole-doped, bottom) BaFe2As2 compounds, (a, d, g) Surface plots of the real (left) and 
imaginary (right) parts of the Lindhard susceptibility within the L = and L = 1 planes, (b, e, h) Respective profiles of /oCQ? ^) along the 
high-symmetry directions, plotted at L = 0, 1/2, and 1. (c, f, i) L -dependence of /o(Qj ^) along the (||^) symmetry axis and at the 
incommensurate peak positions (for doped compounds only). 



with the respective profiles along high-symmetry directions. 
Already in the parent compound, despite the commensura- 
bility of the nesting, a significant in-plane anisotropy of the 
AFM peak is observed both in the real and imaginary parts 
of Xq^ preserving its transverse elongation at all L. This 
clearly indicates that the 42 /m screw symmetry is not to be 
expected in the spin-fluctuation spectrum. In other words, 
our calculations are consistent with the lowered symmetry 
of the spin response that corresponds to the unfolded BZ 



of the Fe-sublattice, as sketched in Fig. 2 (b) . It should be 
emphasized that the asymmetry of the calculated Lindhard 
function along the XT and XZ lines appears only if the ma- 
trix elements of the perturbation are properly taken into 
account in Eq. 1. If the matrix elements are neglected, 2'o(Q) 
becomes four-fold symmetric with respect to the rotation 
around the (||L) axis. 

The stronger response along the transverse direction is 
present at all L values, resulting in an almost vanishing 
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Fig. 5 (color online). Lindhard function (left) and renormalized RPA spin susceptibility (right) in the static limit (co ^ 0), calculated from 
a 3D tight-binding model.^^^ (a) Lindhard function Xo(H,K,0) and Xo(H,K, 1) for the parent (undoped) compound; (b) Corresponding 
renormalized RPA spin susceptibilities jRpA(H,iC, 0) and /Rp^(H,iC, 1) calculated for (7 = 0.8 and J = 0.25 U. (c, d) The same for 7.5% 
electron-doped compound within rigid-band approximation. 



L -dependence, except for the weak intensity modulation 
that is best seen in Fig. 4(c). Due to this modulation, 
Re;^o(Q^ co) — the function that is responsible for the SDW 
instability — is ~1.4% larger at L = 1 than at L = in 
undoped BaFe2As2, which is sufficient to explain the out- 
of-plane component of the 3D AFM ordering wavevector 
(||l) that otherwise cannot be understood using simple 
geometrical nesting considerations. 

As the system is doped by either electrons [Fig. 4 (d-f) ] 
or holes [Fig. 4 (g-i)], the in-plane anisotropy of the (||l) 
peak and, consequently, the absence of the 42 /m symme- 
try become even more apparent. The nesting peaks in the 
Lindhard function develop an incommensurability along 
the directions transverse or longitudinal to Q, respectively, 
which becomes well resolved only at sufficiently high doping 
levels. In the Co-doped compounds below or at the optimal 
doping, where most of the available INS experiments were 
performed, the incommensurability only leads to an addi- 
tional broadening of the peak in the transverse direction, 
and to an increase in the anisotropy ratio as compared to 
the undoped compound. 



The L -dependence of Re;^o ^^ ^^e wave vector i^L) cor- 
responding to stripe-like AFM correlations in the ab plane 
is strongly affected by doping. In undoped BaFe2As2, the 
maximum of Re;^o is found close to 1=1, indicating that 
AFM correlations between Fe layers are favorable [Fig. 4 (c)] . 
Electron doping suppresses the variation of the susceptibil- 
ity along the (||L) line. Figure 4 (f) shows, however, that 
the L -dependence at the maximum of Re;^o^ i-^-^ along the 
(0.56, 0.44, L) line, becomes more pronounced. Hole doping 
[Fig. 4 (i)] leads to even stronger suppression of spin cor- 
relations with Qafm so that Re;^o ^^(22-'^) becomes lower 
than at (||0). The L dependence at the maximum of Re;^o 
at Q = (0.41,0.41,1) is negligible, and only at the local 
maximum Q = (0.56,0.44,1), AFM correlations between 
the layers are still preferable. 



B. RPA spin susceptibility from a 3D tight-binding model 

In order to go beyond the bare spin susceptibility and 
account for electronic interactions, we apply the random 



phase approximation (RPA) to the 3D tight-binding (TB) 
model introduced in Ref. 116, which effectively parameter- 
izes the unfolded DFT band structure calculated for the 
experimental atomic positions. ^^^ Here the Lindhard func- 
tion is calculated from the multiorbital susceptibility^ ^^'^^^ 



(rf,'(Q,«>) = 



-y, 



a;(k)aP*(k)ag(k+Q)<*(k+Q) 
x[/(£,(k+Q))-/(£/k))], (2) 



where p, q, s and t are orbital indices, /i and v label the 
energy dispersion Ey(k), and /(£) is the Fermi function. 
With the summation over all momenta in the first BZ, the 
full 3D dispersion is taken into account. The underlying 
symmetry of the crystal (including the orbital composition of 
the bands) is reflected both in the TB band dispersions E^ (k) 
and in the matrix elements cz^(k), connecting the band and 

orbital spaces. ^^"^ Since there are indications that electronic 
correlations in the iron arsenide systems are moderate, as 
compared to the high-T^ cuprates/^^'^^^ we have included 
the Coulomb repulsion U and the exchange splitting J on the 
Fe sites in the framework of the RPA. Here the multiorbital 
susceptibility of the interacting system is given by^^'^ 

ixr):: = ^x^v + ^xtz ^k^j:. ^xox:, (3) 

where U^^^^ is the interaction matrix in orbital space as 
defined in Ref. 116. In Fig. 5, the Lindhard function 



Xo(Q,a>)=-2,ao)rr(Q,^) 

S = t 

p=q 

and the total RPA spin susceptibility 



(4) 



;^rpa(Q. co) = l E^^D- (Q^ ^)^ (5) 



s=t 
p=q 



calculated for U = 0.8 and J = 0.25 U, are shown in the 
static limit within Q = (Hi^O) and Q = (HKl) planes both 
for the electron-compensated parent compound and for 
the 7.5% electron doping that results from a rigid-band 
shift of the TB bands by 33.5 meV The Lindhard functions 
presented here are not strictly equivalent to those in Fig. 4, 
as they are derived from independent DFT band structures 
and are calculated from a TB fit to the unfolded electronic 
bands, whereas those in Fig. 4 originate directly from DFT 
calculations performed in the backfolded (bet) unit cell. This 
results in subtle differences, such as a sharper nesting peak 
in Fig. 4, that are not essential for the purpose of the present 
paper. We also note that in contrast to Ref. 116, we have 
determined the doping level from the electron count within 
the tight-binding model to ensure internal consistency. The 
notation in Fig. 5 corresponds to the backfolded tetragonal 
BZ and therefore also differs from that of Ref. 116. The 
RPA approach allows for a qualitative analysis of the Q- 
dependence of the measured susceptibility and correctly 
reproduces the location of the signal in the phase space 
and its anisotropy For a quantitative comparison, which 



is outside the scope of the present paper, approximations 
going beyond a standard RPA with momentum-independent 
interactions might be necessary. 

At both doping levels, the dominant feature in Xyipa is 
located around the Qafm wavevector, originating from the 
nesting of hole- and electronlike FS sheets. Its maximum 
appears at a nearly commensurate position in the parent 
compound, but the incommensurability increases drastically 
upon doping as a natural consequence of the rigid-band 
approximation. This is at variance with experiments that 
found a commensurate spin response in a wide range of 
electron doping levels. "^^'^^ This lack of correspondence indi- 
cates that the rigid-band approximation cannot fully account 
for the doping effects in iron arsenides, as suggested earlier 
in several theoretical works. "^^'^^^"^^"^ 

On the other hand, the symmetry of the magnetic spec- 
trum, as well as the tendency to larger anisotropy with 
increased doping, are well captured by the TB model. The 
Lindhard function shows good qualitative agreement with 
the directly calculated one from section II A. The suscep- 
tibility patterns are incommensurate along the transverse 
direction both at L = and L = 1, and therefore do not 
possess the 42 /m symmetry. The RPA renormalization con- 
siderably enhances Im;^o(Q^ co)/co around the nesting vec- 
tor, whereas the strong peak at the F point is considerably 
suppressed due to a much smaller Stoner factor. As a result, 
the overall agreement with experimental spectra that consist 
of a single pronounced feature centered at (||l) is further 
improved. 

In summary, the results of our theoretical calculations 
indicate that the normal-state spin susceptibility contains 
all essential ingredients that are necessary to understand 
the symmetry of the measured INS spectra, both in the 
normal and SC states, on a qualitative level. These include 
both the out-of-plane modulation of the Lindhard function, 
peaked at the Qafm wavevector, and the in-plane anisotropy 
of the nesting-driven peak, which preserves its transverse 
elongation at all L values. Both effects lead to the absence of 
the 42/^ screw symmetry in the spin-excitation spectrum. 

III. EXPERIMENTAL RESULTS 

A. Sample description and experimental details 

The single crystals of BaFe^ 9iNio.o9As2 (T^ = 18K, m ?^ 
4g) and BaFe^ 85C00.15AS2 (T^ = 25K,m^lg) were grown 
by the FeAs-flux method^^"^'^^^ and characterized by energy- 
dispersive x-ray analysis, SQUID magnetometry and single- 
crystal neutron diffraction using the E2 flat-cone diffrac- 
tometer at the Helmholtz-Zentrum Berlin fiir Materialien 
und Energie. Magnetization measurements on several small 
pieces of each sample revealed sharp SC transitions at 
T^ = 18 K and 25 K, respectively, as shown in Fig. 6(a). 
Both in the (HHL) and (HKO) planes, neutron diffraction 
patterns exhibit well defined Bragg spots with narrow mo- 
saicity < 1° [Fig. 6 (b)] and no signatures of multiple single- 
crystalline grains, but with some polycrystalline contamina- 
tion originating both from the main phase and to a lesser 
extent from traces of the (Fe,Co)As flux [see Fig. 6 (c)]. We 
therefore had to optimize the scattering conditions in our 
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Fig. 6 (color online) . Characterization of the samples used for the 
present study, (a) Magnetization curves measured in the magnetic 
field of 10 Oe, applied in-plane, after cooling in the field (FC) and 
in zero field (ZFC). Insets show photos of the samples, (b) Rocking 
curves measured on the (004) reflection in the (HHL) scattering 
plane with a triple-axis spectrometer, (c) Neutron diffraction pat- 
tern of the BaFe^ 9iNio.o9As2 sample in the (HHL) scattering plane. 
Powder lines coming from the Al sample holder and traces of the 
(Fe,Co)As flux are marked by the arrows. 



(a* = b* = 2n/a, c* = 2n/c) that would correspond to a 
simple tetragonal unit cell with the same dimensions. The 
room-temperature lattice constants are a = b = 3.94 A, 
c = 12.86 A for BaFei.9iNio.o9As2 and a = b = 3.92 A, c = 
12.84 A for BaFci 85C00.15AS2. For the sake of a compact 
notation we will seth= 1 in the following and quote the 
energy transfer co in meV 



B. L -dependence of the INS spectra 

In Fig. 7, several representative longitudinal Q-scans 
across the AFM wavevector are shown. "^^ One can see that 
both in the normal and SC states, the signal is well fitted by 
a single Gaussian peak with a linear background, showing 
no signatures of incommensurability along this reciprocal- 
space direction within the low-energy range of up to ~ 2 A. 
In Fig. 8, we show the energy dependence of the experi- 
mentally measured imaginary part of the spin susceptibility 
X^XQ, co) at Q = (||l) for both samples at even and odd L, 
obtained from the raw INS data after background subtrac- 
tion and Bose-factor correction. The measured signal has 
also been corrected to account for the energy-dependent 
fraction of higher-order neutrons. The data were acquired 
by performing a series of full Q-scans similar to those shown 
in Fig. 7 at different fixed energies and an energy scan at 
Q = (||l). To estimate the background for the latter, we 
used a linear interpolation for the background obtained 
from Gaussian fits to the full Q-scans, or measured points 
appropriately offset to both sides from i^^l^L). The error bars 
correspond to one standard deviation of the neutron count 
and do not include the normalization errors. The two left 
panels of Fig. 8 show data on BaFe^ 9iNio.o9As2, measured in 
the SC and normal states at L = 1 and 3 [panel (a)] and at 
L = 2 [panel (c)]. The respective data for BaFe^ 85C00.15AS2 
are shown at the right. 



INS measurements by avoiding the appearance of spurious 
inelastic peaks caused by such contamination. No structural 
or SDW transitions were detected down to 2 K in both sam- 
ples, consistent with the known phase diagrams. "^^'^^^"^^^ 

The INS measurements were performed at the triple-axis 
spectrometers PANDA and PUMA (FRM-II, Garching), IN8 
(ILL, Grenoble), and 2T (LLB, Saclay). The instruments 
were operated in their high-flux setup without collimators, 
using focussed pyrolytic-graphite (002) monochromators 
and analyzers. Measurements were done in the constant-Zcf 
mode, with k^ = 1.55 A"^ (Ef = 4.98 meV) or kf = 2.662 A"^ 
(Ef = 14.7 meV). Correspondingly, either a cold Be-filter or 
two pyrolytic-graphite filters were used for higher-order 
neutron elimination. 

The data for the present work were collected in the (HHL) 
and (^HK [H-\-K]) scattering planes. Throughout this paper 
we are using backfolded tetragonal notation, ^^^ in which 
Qafm — (2 2 -'■) corresponds to the AFM ordering wavevector 
of the parent compound. We quote the wavevector Q = 
(^HKL) in reciprocal lattice units (r. l.u.), i.e. in units of 
the conventional reciprocal lattice vectors a*, b*, and c* 
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Fig. 7 (color online). Several raw Q-scans for BaFei85Coo.i5As2, 
measured along the longitudinal direction in the SC state (top row, 
r = 4 K) and in the normal state (bottom row, T = 60 K) at three 
different energies: Smey 9.5 mey and 16meV The solid lines 
represent Gaussian fits with a Hnear background. The background 
is indicated by dashed lines. 
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Fig. 8 (color online) . Imaginary part of the spin susceptibility 
at odd (top) and even (bottom) L in the normal and SC states. 
The left column shows data for BaFei9iNio.o9As2 at Q = (||l) 
and (||3) in (a) and at (||2) in (c). The right column shows 
corresponding data for BaFe^ 85C00.15AS2. The data points were 
obtained from constant- co scans and constant-Q scans, as described 
in the text. The solid lines are guides to the eye. Different symbol 
shapes represent data obtained in different measurements. 



Already in the normal state, a difference between odd 
and even L values can be observed. For both samples, the 
normal-state spectral weight, integrated over Q and co up to 
14mey is ~ 60% larger at odd than at even L. Such a differ- 
ence cannot be a consequence of the magnetic form factor, 
which would be smaller at L = 1 than at L = 0, producing 
the opposite effect. On the other hand, this difference is rem- 
iniscent of the SDW phase of the parent compounds, where 
low-energy magnon branches are present only near mag- 
netic Bragg peaks at odd L, whereas spin waves at even L 
are gapped and thus yield zero intensity at low energy.^^"^^ 
However, in the paramagnetic state, the normal-state inten- 
sity at even L is only moderately suppressed [cf. Fig. 4 (c) 



and (f)]. Here we note that the absence of any magnetic 
Bragg intensity at (| - 1 0) or (-|| 0) in the SDW state is 
fully consistent with the unfolded-BZ scheme. Indeed, as 
can be seen from Fig. 3, these two X points correspond to 
the zone center in the doubly-folded magnetic BZ, which 
means that the influence of the As superstructure would lead 
to an appearance of magnetic Bragg-peak replicas at these 
points. In a twinned crystal, this would imply equivalency 
of all the (±1 zb I L) points up to the magnetic structure 
factor. The fact that these replicas have not been observed 
by neutron diffraction indicates that the structure factor for 
the As-superstructure reflections is negligibly small or zero. 
In other words, no folding of the magnetic signal occurs 
due to the As sublattice, and hence the unfolded-BZ scheme 
is perfectly justified. Our results presented in this and the 
following sections serve to generalize these arguments to 
the inelastic magnetic signal. 

At first, we consider the low-temperature spectra that 
exhibit the SC resonant mode. We define the resonance en- 
ergy cOj-es ^s the maximum of x^XQ^ ^) ii^ the SC state and 
discriminate between its value at even and odd L, C0res,even 
and cOj-es^odd^ where necessary. The dashed vertical lines 
mark these positions for odd and even L in the upper and 
lower panels of Fig. 8, respectively. We note that the res- 
onance energies at odd and even L differ by more than 
2meV in BaFei9iNio.o9As2, in agreement with Ref. 86. In 
contrast, this difference is only about 1 meV in optimally 
doped BaFci 85C00.15AS2, as seen in Fig. 8 (c) and (d). 

Fig. 9 shows the magnetic intensity evolution near the 
resonance energy along L, obtained from Gaussian fits of 
full constant-energy scans around Q = (||^) ^t 8 meV 
Similarly to the normal-state intensity, it is modulated pe- 
riodically in L (up to the magnetic form factor), analogous 
to the magnons in the parent compound. Two factors can 
be responsible for the observed modulation. First, as the 
normal-state intensity is already lower at even L, it will 
preserve this modulation after redistribution of the spectral 
weight due to the opening of the SC gap below T^. Second, 
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Fig. 9. L -dependent magnetic intensity of BaFei9iNio.o9As2 in the 
SC state atQ = (||L) and 8 meV (close to the resonance energy). 
The dashed line shows the Fe^^ spin-only magnetic form factor. 
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Fig. 10. Comparison of momentum profiles at even and odd 
L at fixed energies that are below co^^ for even L, but above 
it for odd L. (a) BaFei9iNio.o9As2, T = 3K and co = 3meV 
(b) BaFci 85C00.15AS2, r = 4K and CO = 4meV 
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the higher energy of the resonance at even L is closer to (or 
even within) the particle-hole continuum, which may result 
in stronger damping and additional intensity reduction. 

Not only cOj-es and the spectral weight of the resonance, 
but also the energy range below the resonance peak that 
is depleted upon entering the SC state (which we refer to 
as the SC spin gap) depends on L. We define the spin- 
gap energy co^g as the intersection of the low-energy linear 
extrapolation of x^XQ, ^) ^t 2 or 4 K with the x^^ = line 
(Fig. 8). Inspection of Fig. 10, where we compare constant- 
energy scans at even and odd L for both samples, clearly 
shows that co^^ is larger at even L . We remark that the SC 
spin gap should not be mistaken for the SC gap A to which 
it is only indirectly related: co^g is determined by the energy, 
cOj-es. and the width of the resonant mode. 

Recalling that the X points in the BZ for odd and even 
L values are equivalent due to the above-mentioned screw 
symmetry, we can now conclude that this symmetry is ab- 
sent in the spin-excitation spectra of both samples based 
on the following evidence observed in the out-of-plane 
direction: (i) different normal-state intensities at odd 
and even L; (ii) different resonance energies cOres,odd ^^^ 
cOj-es^even^ (iii) periodic L -dependent intensity of the reso- 
nance; (iv) the corresponding difference of the spin gaps 

^sg,odd ^^' ^sg,even* 

C. In-plane anisotropy of the spin susceptibility 

Another piece of evidence for the lowered reciprocal- 
space symmetry is associated with the in-plane anisotropy 
of the measured INS intensity.^^"^^'^^'^^ In Fig. 11 (a), we 
show experimental constant-energy maps, interpolated from 
a series of triple-axis Q-scans in the vicinity of (||l) and 
(— 1|0) wavevectors, measured in the (HK [H-\-K]) scatter- 
ing plane. We compare them with the calculated dynamic 
spin susceptibilities of the paramagnetic tetragonal phase, 
plotted in the equivalent regions of Q-space surrounding 
the X points. Panel (c) shows the imaginary part of the 
Lindhard function lmxoiH,K,0) (left) and lmxoiH,K,l) 
(right) in the vicinity of Qy and Qafm? respectively, for 7.5% 
Co-substitution, as calculated by DFT in the virtual crystal 
approximation. Panel (d) displays the respective results for 
the RPA-enhanced susceptibility Im^'RPAC'^ + Qx? '^ + Qy^ 0) 
and Im;^RpA('Ti + q^, u + q.^, 1) [same as in Fig. 5 (d)], calcu- 
lated in the rigid-band approximation from the TB model^^^ 
at 7.5% electron doping. 

Notably, the transverse elongation of the susceptibility 
pattern is preserved at all L values both in the measured 
INS signal and in the results of both calculations, meaning 

that the longer axis of the ellipse is oriented either along XZ 

or along Xr directions for even and odd L, respectively. This 
anisotropy is insensitive to the SC transition and persists 
also in the normal state. Neither the widths of the peaks nor 
their anisotropy experience any change across T^ within our 
experimental accuracy, as evidenced by Fig. 1 1 (b) . 

In comparison to the magnetically ordered parent 
compound, ^"^^ which exhibits a steep spin wave disper- 
sion cone around Qafm^ ^s shown by a small black ellipse 
in Fig. 11 (a), electron doping tends to increase the trans- 
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Fig. 11 (color online), (a) Experimental intensity distributions 
for BaFe^ 85C00J5AS2 near Qy (left) and Qafm (right), measured in 
the (HK[H -\-K]) scattering plane in the SC state (T = 4K) at 
the resonance energy (9.5 meV). The small black ellipse around 
( II 1) is a 9.5 meV cross-section of the spin wave dispersion for the 
CaFe2As2 parent compound,^^'^^ shown for comparison. The white 
dotted lines are BZ boundaries, (b) Comparison of the longitudinal 
(LO) and transverse (TR) cross-sections of the data from panel (a) 
around L = (left) and L = 1 (right) . (c) The Lindhard function 
ImjDFT(^)/^ ^t 7.5% Co-doping, calculated by DFT in the same 
reciprocal space regions, (d) The same for the RPA-renormalized 
low-energy spin susceptibility Im/RpA(ct))/(x) [same as in Fig. 5 (d)], 
calculated from a TB model in the rigid-band approximation. 
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Fig. 12 (color online), (a) Longitudinal (LO) and transverse 
(TR) widths of the commensurate peaks around Qafm (^ = 1) 
and Q|| (L = 0) at 9.5 meV versus temperature, (b) The same 
widths versus energy transfer at low temperatures. Solid lines 
are results of a global fit to the data in both panels (see text) 
using Eq. (6). Resolution-corrected dependencies are shown 
by dashed lines, (c) The resolution-corrected anisotropy ratio 
A = (w^R — Wlo)/(wtr + Wlo) as a function of temperature, com- 
pared to the respective values for the magnetically ordered^^'^^ 
and paramagnetic^^ states of CaFe2As2. (d) The same ratio as a 
function of energy transfer. The dashed line gives the anisotropy 
from the global fit. The dotted line is derived from the high-energy 
dispersion reported for BaFci 87C00.13AS2 in Ref. 30. 



verse incommensurability [cf. Fig. 4 (c, d)] and, in addition, 
leads to softening of spin excitations predominantly in the 
transverse direction.^^ This results in a rapid increase of 
the anisotropy ratio w^ith increasing doping. The emerg- 
ing pattern resembles the "unusual quasi-propagating ex- 
citations" observed at higher energies in a similar com- 
pound by Li et al,^^ as v\^ell as the pair of incommensurate 
peaks seen in FeTei_^Se^ (Refs. 34,91,92,94,95). In the 
light of our present results, the former can be understood 
as tw^o incommensurate branches of itinerant Stoner-like 
excitations, driven by FS nesting, as in the case of iron 
chalcogenides."^'^'^^'^'^'^^ The fact that such incommensura- 
bility has not been resolved experimentally at low^ energies 
is not surprising, because for sufficiently small doping levels 
at w^hich the overwhelming majority of INS experiments 
w^as performed, the two incommensurate peaks merge into 
one due to their finite width, resulting in a broad commen- 
surate peak elongated in the transverse direction. Similar 
measurements of strongly overdoped samples are therefore 
necessary to confirm this scenario and the emerging similar- 
ity to the 11 -compounds. 

In order to quantify the observed in-plane anisotropy and 
compare it with previous experiments, in Fig. 12 we plot 
the temperature and energy dependence of the measured 
full width at half maximum (FWHM) of the commensurate 
inelastic peak along the longitudinal (LO) and transverse 
(TR) directions for L = and L = 1. In the longitudinal 
direction, the resolution-corrected width of the peaks Wlq 
(dashed line) was already quantified for the same sample 
by a fit to the Moriya formula^"^^ in Ref. 76. To extract the 



anisotropy ratio A = (w^r — Wlo)/(wtr + v^^o), we have 
fitted the experimentally measured FWHM of the peaks 
in the longitudinal (Wlq) and transverse (W^r) directions 
(solid lines in Fig. 12) using the following equations: 



Wlo(co, r) = ^wl^{co,T)+R^', 



Wtr(co, T) = 



1+A 
1-A 



Wlo(co, T) 



(6) 



+ R\ 



The fitted value of the effective resolution, R = 0.066 =b 
0.004 r. I.U., was used to perform resolution correction of 
the experimental data and calculate the anisotropy ratio 
that is presented in panels (c) and (d). By setting the reso- 
lution to a constant, we relied on the fact that the calcu- 
lated instrumental resolution is nearly isotropic and does 
not vary within our region of interest by more than ~ 10%. 
The effective momentum-space resolution resulting from 
our fit (hatched region in Fig. 12) is somewhat lower than 
the calculated instrumental resolution (i^min ^ 0.04 r. l.u.). 
The difference may indicate a finite-size limit on the fluc- 
tuating domains imposed by the random distribution of 
dopant atoms and/or a slight inhomogeneous broadening 
due to variations of the doping level across the sample. With 
this reasonable assumption, the entire data set can be de- 
scribed by a single, temperature- and energy-independent 
anisotropy parameter. A similar fit based on the instrumen- 
tal resolution alone (without finite-size or inhomogeneous 
broadening) would yield an anisotropy parameter that in- 
creases with temperature, which would be highly unusual. 
The anisotropy ratio A = 0.41 ± 0.02 that results from 
the global fit to our data is shown in Fig. 12 (c) and (d) by 
the dashed line. This value corresponds to the aspect ratio 
w^r/wlo = 2.4 =t 0.1, which is nearly a factor of 2 larger 
than the respective ratio of spin wave velocities (~1.4) 
in the undoped CaFe2As2, according to Refs. 61 and 62. 
The dotted line in Fig. 12 (c) shows that the anisotropy 
ratio remains nearly constant across the SDW transition, 
as estimated from the paramagnetic-state data measured 
at r = 180 K by Diallo et al^^ On the other hand, the 
anisotropy ratio of 0.44 extracted from the high-energy TOF 
data on a similarly doped BaFe^ 87C00.13AS2 compound"^^ 
[dotted line in Fig. 12 (d)] perfectly coincides with our value. 
This agreement confirms the energy independence of the 
anisotropy and indicates that the difference in the peak 
widths originates mainly from two unresolved incommensu- 
rate peaks, in agreement with our DFT calculations, rather 
than from an anisotropic broadening caused by the finite 
correlation lengths of the spin excitations. ^^'"^^ Despite the 
present lack of Q-resolved INS data on hole-doped com- 
pounds, the results of our susceptibility calculations from 
section II A allow us to predict that the anisotropy of the 
spin-excitation spectrum should vanish and subsequently 
switch to the longitudinal orientation as the system is doped 
with more holes. 

D. Doping dependence of the resonance 

In order to investigate the doping dependence of the res- 
onance and its L -modulation, we summarize in Fig. 13 (a) 
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Fig. 13 (color online), (a) Doping dependence of co^^^ at odd and even L in BFNA and BFCA studied here (full symbols) and in previous 
works (empty symbols), as referred to in (a) and (c). The blue line follows the average T^, rescaled to 4.3 at its optimum. ^^^"^^^ The red 
line is a guide to the eye. (b) Linear extrapolation of the resonance intensities to the energy axis, as compared to the onset of particle-hole 
continuum. The hatched region covers the range of directly measured 2A values for the larger gap in nearly optimally doped BFCA, 
estimated by various experimental techniques. ^^^"^^^ (c) Resonance energy versus T^. ^p^ for different Fe-based superconductors. For the 
compounds with dispersing resonance, only co^^^ at odd L is shown. 



our results together w^ith other studies of electron-doped 
BaFe2As2 (no momentum-resolved data for hole doping^^^ 
are available so far).^^"^"^ To put the cOj-es values from dif- 
ferent compounds and doping levels on the same scale, we 
divided cOj-es by the optimal k^ T^^opt ^^^ normalized the dop- 
ing level by the optimal doping level, respectively. While 
<^res,odd valucs (bluc symbols) inQ = (||L) fall onto the 
blue dotted line which follows the average T^ in the phase 
diagrams from Refs. 126-129, cOj-es^even values (red symbols) 
do not follow Tc, but rather stay at higher energies than 
^res,odd ^^ ^^c uudcrdopcd region, in agreement with a 
similar recent study. ^^^ As a consequence, the difference 
between cOres,odd ^^^ ^res,even incrcascs with underdoping 
(as can also be seen in Fig. 8). 

The integrated intensity of the resonance is influenced, 
in particular, by its proximity to the particle -hole contin- 
uum with an onset at 2 A. As a consistency check, we 
therefore plot in Fig. 13 (b) the Q- and co -integrated inten- 
sities of the resonance at odd and even L versus its energy. 
Since in an RPA description the spectral weight of the reso- 
nant mode is roughly proportional to its excitonic binding 
energy, ^^^'^^^'^^^ under the assumption of L -independent 
onset of the particle-hole continuum, a linear extrapolation 
of the two intensities onto the energy axis gives us a rough 
lower estimate of 2 A — the point where the resonance inten- 
sity is fully suppressed by particle-hole scattering (for similar 
analysis in cuprates, see Ref. 102). For the Co-doped com- 
pound, such an extrapolation results in 2Abpca ^ ll.Smey 
which indeed falls in the middle of the range of values re- 
ported from direct measurements^^^"^"^"^ (hatched region). 
Since SC gap measurements for the Ni-doped compound 
are scarce, we resort to calculating the coupling constant 
2A/k^T^ = 6.8 that results from the extrapolated gap of 
2Abpna ^ 10.6 meV On the one hand, it agrees with the uni- 
versal value of 7 =b 2 that was reported for the larger gap in 
various two-gap ferropnictides^^^ and coincides with that of 
6.8 (or 6.6) derived from combined ARPES and muon-spin 



rotation (/iSR) 



144,145 



and specific-heat^^ measurements on 



Bai_^K^Fe2As2, respectively. On the other hand, it exceeds 
the maximum coupling constant of 2A/A:Brc ^ 5.0 that was 
recently inferred^"^^'^"^"^ from specific-heat measurements on 
BFCA. The non-linear dependence of the larger gap on T^, 
reported in Ref. 137, would result in a much lower esti- 
mate for 2Abpna ^ 6-9 meV in the Ni-doped sample, under 
the assumption that this dependence is universal among 
122-compounds. Such low value would imply a consider- 
able overlap of the resonance peak with the particle-hole 
continuum, which could explain its broad width in energy. 

The successful application of the simple scaling relation 
with L -independent particle-hole continuum indicates that 
the distance between the resonance and the continuum 
2 A — cOj-es is L -dependent, as otherwise the agreement with 
directly measured gap values would be coincidental. In 
other words, the L -dependence of the resonance energy and 
intensity alone does not necessarily imply a A:^ -dependent 
energy gap, as suggested previously,^^ but more likely is 
a natural consequence of the normal-state intensity mod- 
ulation. While a SC order parameter that differs at odd 
and even L values is conceivable and was even supported 
by experimental evidence, ^^'^^'^"^"^ it can only result from 
the normal-state properties of the "pairing glue", and thus 
does not appear to be the primary reason for the dispersing 
resonant mode. 

Finally, in Fig. 13 (c) we combine our data with all the 
previously reported data^'^"^^ to show how the resonance 
energy cOj-es ^^ odd L depends on T^. Filled symbols are 
extracted from the present work, whereas empty symbols 
are from the references indicated in the plot. We see that 
the resonance energy scales linearly with T^ from the lowest 
to the highest critical temperatures with a universal ratio 
of cOi.es,odd/^B^c ^ 4.3 for all hole- and electron-doped Fe- 
based superconductors investigated so far. It is interesting 
that such simple linear scaling apparently breaks down at 
even L values, ^^^ as also seen in Fig. 13 (a). 
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IV DISCUSSION AND CONCLUSIONS 

A. Normal state 

As we have demonstrated, the eUiptical shape of the spin 
excitations within the L = const planes shows no measur- 
able L -dependence (apart from an intensity modulation) 
and is insensitive to the SC transition. Therefore, the origins 
of this anisotropy are to be found in the properties of the 
normal (paramagnetic) state. An anisotropic spin correla- 
tion length that is larger in the direction parallel to the AFM 
propagation vector than in the transverse (ferromagnetic) 
direction has been proposed as a possible explanation. ^^^'^^^ 
Although such description is successful in the low-energy 
region, where the two spin wave branches are not resolved, 
it clearly fails to describe the anisotropic spin wave veloc- 
ities that become evident at higher energies in the para- 
magnetic state,"^^'"^^ mimicking the behavior of the parent 
compounds. ^^ This implies that the larger momentum width 
of the spectrum in the transverse direction is more likely 
to be a result either of two unresolved spin-wave branches 
that are less steep than the longitudinal ones, or of the 
incommensurability of the nesting peaks at co = 0. The 
results of our DFT calculations support the incommensu- 
rate nesting scenario, similar to that inferred earlier from 
nuclear-magnetic-resonance measurements ^"^^ and to the 
one proposed for the iron chalcogenides.'^'^'^^'^^'^^ In such a 
case, the anisotropy results from FS nesting, and not from 
an "electronic liquid-crystal state" that arises spontaneously 
from electron-electron interactions. ^"^'^^""^^'^"^ The latter state 
has been invoked for the cuprates based in part on the strong 
temperature dependence of the in-plane anisotropy of the 
spin excitations,^^ which is not observed in the 122-iron- 
arsenide system (Fig. 12). If our prediction of the rotated 
(longitudinally elongated) susceptibility profile in the hole- 
doped compounds were confirmed experimentally, it would 
provide additional support for this scenario. 

Although the exact causes for the anisotropy are still un- 
der debate, it has been argued that remnant magnetism 
persisting above the AFM transition, in the form of fluctuat- 
ing magnetic domains^^ or the above-mentioned "electronic 
nematic" ground state, is involved. However, as we have 
demonstrated, the spin-fluctuation spectrum possesses con- 
siderable anisotropy and does not fully follow the crystallo- 
graphic symmetry even in the normal (paramagnetic and 
tetragonal) state, when no electronic nematicity is assumed. 
Therefore, the observed in-plane anisotropy in the doped 
compounds does not necessarily imply a symmetry-broken 
ground state, but has a more trivial structural origin. As 
the primitive structural unit cell of the 122-compounds con- 
tains two Fe atoms, its size is twice larger, as compared to 
that of the Fe-sublattice. Because the magnetic INS signal 
originates predominantly from the latter, with no magnetic 
moment being induced on the As sites,^"^'^^ the symmetry of 
the spectrum is determined by the unfolded BZ. In the real 
bet BZ, both the electronic bands and the spin susceptibility 
are folded, but the matrix elements that are responsible for 
the intensities of the primary features and their replica (an 
analog of the dynamic structure factors) are such that no 
abrupt change in the magnetic spectrum can be seen as long 



as the folding potential remains sufficiently weak. Similar 
effects were observed and discussed in relationship to the 
single-particle spectral function (see, for example, Ref. 149 
and references therein), where matrix elements are also 
famous for shaping some of the spectral features. ^^^'^^^ 

The periodic modulation of the magnetic spectral weight 
with L can also be explained by the L -dependence of the 
spin susceptibility observed in our normal-state DFT calcu- 
lations. Although the variation of the Lindhard function 
between L = and L = 1 is weak in the parent compound, 
it can possibly be enhanced by the Stoner-like renormaliza- 
tion effects to an amplitude comparable with experimental 
observations. The maximum of Re;^o(Q^O) in the parent 
compound occurs at Qafm — (| I^) ^^^ ^^^^ determines the 
AFM ordering wavevector. 

In the SDW state, excitations at even L correspond to 
zone-boundary magnons which, due to a combination of 
intra- and interlayer coupling parameters, have a substantial 
gap of ~ 80 meV.^^~^^ In contrast, at high doping levels the 
magnetic response is virtually L -independent.^^ Two mecha- 
nisms are likely to provide the connection between these two 
limiting cases: First, when approaching the magnetically 
ordered state from higher doping levels, the paramagnon 
mode softens at Qafm? ^^^ the in-plane magnetic correlation 
length increases. ^^^ As a consequence, one can expect the 
out-of-plane magnetic correlations to become more efficient 
in stabilizing the mode and its gapped response at even L. 
Second, when starting from the ordered state, the increasing 
damping of the mode with doping will progressively redis- 
tribute spectral weight towards lower energies, including 
the gapped region around even L}^^ At our intermediate 
doping levels we thus observe a moderate L -modulation in 
the normal state [Figs. 8 and 13 (a)]. 

Starting from the paramagnetic state, one sees that the 
similarity of excitation spectra in the magnetically ordered 
and normal states does not imply that the anisotropy of the 
SDW state survives above the structural transition in the 
form of "spin nematic correlations". On the contrary, the 
symmetry-breaking L -modulation is present in the tetrag- 
onal phase for reasons not related to magnetic ordering, 
whereas the SDW instability that occurs on top of the para- 
magnetic state upon cooling or decreasing the doping is 
predetermined by this modulation, so that the AFM prop- 
agation vector coincides with the strongest nesting-driven 
peak in Re ;^o(Q^ 0)- ^^ electronic nematic state also appears 
implausible in view of the temperature independence of the 
in-plane anisotropy (Fig. 12), which is in sharp contrast to 
the strongly temperature dependent, order-parameter-like 
behavior observed in YBa2Cu305+3, (Ref. 20). Our conclu- 
sions about the nonmagnetic origin of the missing symmetry 
are additionally supported by the following evidence: (i) 
experimentally observed enhancement of the anisotropy in 
the doped compound with respect to a magnetically ordered 
parent, which agrees with the increased transverse incom- 
mensurability seen in the DFT calculations; (ii) temperature- 
independence of the anisotropy even in the parent com- 
pound, including its insensitivity to the presence of static 
AFM order. Independently of its origins, the symmetry of 
the normal-state spin-fluctuation spectrum may have im- 
portant implications for the SC order parameter under the 
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assumption of spin-fluctuation-driven superconductivity. It 
was argued, for example, that the transverse elongation of 
the spin-fluctuation profile stabilizes the s"^ pairing state. ^^^^ 



B. Superconducting state 

Now, let us consider the SC state properties and discuss 
the implications of our results for the SC pairing mecha- 
nism. First, we note that while the conventional unit cell 
contains two FeAs layers, the primitive cell, from which the 
BZ is constructed, contains only one. Thus, in contrast to 
cuprates only one resonant mode is expected; the different 
resonance energies at even and odd L are therefore to be 
attributed to an L -dispersion rather than to a mode splitting, 
in agreement with a recent report. ^^^ This dispersion sig- 
nals the non-negligible 3D character of the electronic band 
structure and its importance for the description of the SC 
state. Indeed, there is compelling evidence for such three- 
dimensionality both from ARPES"^^"^^ and band structure 
calculations. "^^'^^ 

Since the calculations we have presented here are limited 
to the normal state, our discussion of the magnetic spectrum 
in the SC state has to remain on a qualitative level. However, 
in view of the normal-state L -modulation, already a mini- 
mal model like RPA is expected to capture the L -dispersion 
of cOj-es '- By virtue of the resonance condition, namely the 
vanishing denominator in ;^o(Q^ ^)/[l ~^(Q);^o(Q^ ^)]5 the 
modulation is carried over into the SC state. Here both the 
bare susceptibility Xo and the interaction / can depend on L. 
We keep / deliberately general — there is no need to refer, 
for instance, to a t-J model,^^'^^^ whose applicability to the 
iron arsenides is being disputed. 

We now put our considerations into a broader context 
by comparing our results to the resonant phenomena in 
YBa2Cu305+-^ (Fig. 14). The latter consists of nearly in- 
dependent CUO2 bilayers and exhibits manifestly 2D elec- 
tronic structure and SC gap. One observes two distinct, 
non-degenerate resonances due to the difference in both 
the bare susceptibility Xo ^^^ ^he interaction / between 
the even and odd channels, ^^"^"^^^ which can be ultimately 
tracked back to the contrast between the intra- and inter- 
bilayer hopping and interaction terms. On the contrary, in 
our iron-arsenide samples, this contrast vanishes and we ob- 
serve a single resonance, which in addition disperses for the 
reasons described above. Thus, both systems represent dif- 
ferent limiting cases of a more general model with coupled 
bilayers and possibly 3D electronic structure. Fig. 14 (b), 
where we expect two resonant modes which both disperse 
and exhibit an intensity modulation along L, depending on 
the effective coupling. 

The similarity between the doping dependence of the 
out-of-plane dispersion bandwidth in the 122-family of iron 
arsenides [Fig. 13 (a)] and the even-odd resonant-mode 
splitting in bilayer cuprates^^^'^^"^ supports our juxtaposition 
of the two systems: In both cases, the even-odd difference 
increases when moving towards the magnetic quantum crit- 
ical point. Whereas the vanishing difference in Fig. 13 (a) 
around optimal T^ suggests that it is determined by the 
proximity to the magnetic instability, emphasizing the im- 
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Fig. 14 (color online). Illustration of the evolution from an 00- 
layer system like BFCA to a bilayer system like YBCO in terms of 
inter- and intra-bilayer distances and effective interactions. For the 
equidistant limit (left), a single dispersing resonant mode is ob- 
served, whose intensity modulation (shown here by the brightness 
of the curve) is mainly governed by the closeness to the particle- 
hole continuum with an onset at 2 A. The dashed line depicts the 
replica that gains intensity only after the equivalency of the layers 
is broken (middle panel) . For alternating interlayer coupling, the 
resonance splits into odd and even modes, which become non- 
dispersive for the case of YBa2Cu305+^ with nearly independent 
bilayers (right). 



portance of the out-of-plane magnetic coupling in the ar- 
senides, recent measurements suggest a persistent even-odd 
difference even beyond optimal doping level, ^"^^ indicating 
that it rather scales with T^. More detailed experimental 
and theoretical work is necessary to settle this point. 

Next, we address the implications of the linear rela- 
tionship between co^^^ and T^ with cOj-es^odd/^B^c ^ 4.3, 
Fig. 13 (c). First, the lower value of this ratio, as compared 
to that of 5.3 for cuprates, supports the notion of a weaker 
SC pairing in Fe-based superconductors. "^^ Second, the valid- 
ity of the linear relationship for all Fe-based superconductors 
hitherto studied, independent of the doping carrier type and 
over the entire studied doping range, '^^'^'^"^^ suggests that 
models that attribute the resonant mode to an excitonic 
bound state within the SC gap may be more straightfor- 
wardly applicable to Fe-based superconductors than they 
are to the cuprates. Whereas in cuprates deviations from the 
linear relationship accompany the increasingly anomalous 
physical properties at underdoping,^^^'^^"^ the resonance in 
Fe-based superconductors is remarkably insensitive to the 
proximate magnetic state and even coexists with it at very 
low doping. ^^^ 

Finally, we remark that while in a recent report a large 
L -dispersion bandwidth was related to the presence of long- 
range magnetic order or pronounced spin correlations, ^^^ 
here we observe appreciable bandwidth (only ~ 35% less 
than in Ref. 105) in a paramagnetic compound, which we as- 
sociate with the normal-state intensity modulation that can 
be qualitatively reproduced even in the Lindhard function 
calculated for the nonmagnetic ground state. 
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